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Abstract 

We consider the charmonium production in thermahzed hadronic 
medium created in ultrarelativistic heavy ion cohisions at LHC en- 
ergy. The calculations for the secondary J/ijj and production by 
DD annihilation are performed within a kinetic model taking into 
account the space-time evolution of a longitudinally and transversely 
expanding medium. We show that the secondary charmonium pro- 
duction appears almost entirely during the mixed phase and it is very 
sensitive to the charmonium dissociation cross section with co-moving 
hadrons. Within the most likely scenario for the dissociation cross 
section of the J/ij: mesons their regeneration in the hadronic medium 
will be negligible. The secondary production of ip' mesons however, 
due to their large cross section above the threshold, can substantially 
exceed the primary yield. 



1 Secondary charmonium production 

The initial energy density in ultrarelativistic heavy ion collisions at LHC 
energy exceeds by a few order of magnitudes the critical value required for 



1 



quark-gluon plasma formation. Thus, according to Matsui and Satz [1,2], one 
expects the formation of charmonium bound states to be severely suppressed 
due to Debye screening. The initially produced cc pairs in hard parton scat- 
tering, however, due to charm conservation, will survive in the deconfined 
medium until the system reaches the critical temperature where the charm 
quarks hadronize, forming predominately D and D mesons. An appreciable 
fraction of cc pairs and consequently D,D mesons produced in Pb-Pb col- 
lisions at LHC energy can lead to an additional production of charmonium 
bound states due to reactions such as: DD* + D*D + D*D* ^ ip + n and 
D*D* + DD — s> ^ + p as first indicated in . In this work we present a quan- 
titative description of the secondary J/ ip and ip' production due to the above 
processes from the thermal hadronic medium created in Pb-Pb collisions at 
LHC energy. 



2 Thermal production kinetics 

The charmonium production cross section a^jj^^i^ can be related to the 
hadronic absorption of charmonium c^h^Dfj, through the detailed balance 
relation ^ 

^DD^^h = djyl){—^)'^a^f^^jy£), (1) 

where kli, = ~ {^a + mbY][s — {nia — mfe)^]/4s denotes the square of the 
center of mass momentum of the corresponding reaction and h stands for p 
or 77, whereas ip for J/ip or mesons. The degeneracy factor d takes the 
values: dj^tj^* = dj:)*[) = 3/4, djj^j^, = d^^j^ = 1/4 with the pions in the final 
state and c?^^* = dj^*i) = 9/4, dj^i, = 27/4, rf^*^, =3/4, for the processes 
with p mesons. 

The magnitude of charmonium absorption cross section on hadrons is 
still, however, theoretically not well under control. There are four models 
which have been proposed to give an estimate for the J/ ip dissociation cross 
section on pions: (1) the constituent quark model (2) the comover model 
[|, P], (3) an effective hadronic Lagrangian and (4) a short distance QCD 
approach ^j. 

In the framework of the constituent quark model the J/ip absorption 
on a pion is viewed as a quark exchange process where the charm quark 
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Figure 1: The cross section for J/ip + -k DD* + D*D absorption as a 
function of energy: curve (1) is calculated within a non-perturbative quark 
exchange model [4]; curve (2) corresponds to the constant cross section of 3 
mb [5,6]; curve (3) was obtained within the framework of a meson exchange 
model [7]; curve (4) is derived in the context of perturbative QCD [8]. 

changes its side with a light quark. The cross section for the tt + J/ip — »• 
D*D + DD* + D*D* process, in this approach, was calculated in the non- 
relativistic quark model in terms of the first Born approximation The 
energy dependence of the above cross section for the final state D*D + DD* 
is shown in fig. 1. One sees that the cross section abruptly increases just 
above threshold and reaches its peak value of about 6 mb. The large value 
of the cross section is mostly due to the particular modelling of the confining 
interaction between the quarks which is taken as attractive, independent of 
the colour quantum number of the affected quark pair. 

The large charm quark mass, rric ~ 1.2 — 1.8 GeV and J/ip binding energy, 
2Me, — Mj/^ ~640 GeV in comparison to Aqcd, as well as small size of J/ip, 
rj/^ ~ 0.2fm<< Aqqjj, have been used as an argument to calculate the char- 
monium dissociation on light hadrons in terms of a QCD approach P, |^. The 
J/ip — h cross section in this approach is expressed in terms of hadronic gluon 
distribution functions, making use of the short distance QCD method based 
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Figure 2: Model for time evolution of temperature in the expanding QCD 
medium crated in Pb-Pb collisions at LHC energy. 



on sum rules derived from the operator product expansion. The resulting 
energy dependence of the J/ip absorption cross section on pions calculated 
within QCD |^ is shown in fig. 1. The hadronic gluon distribution functions 
are strongly suppressed at high gluon momenta. Thus, since the dissociation 
of J/ip requires hard gluons, the absorption cross section J/ipn — > DD be- 
comes very small just above threshold as seen in fig. 1. Recent analysis of 
J/ip photoproduction data confirms the relation between the charmonium- 



hadron interactions and the hadronic gluon structure p, |T0|, |TT| indicating a 
small value of the J/ absorption cross section on hadrons. 

The scattering of J/ip on tt ot p meson can be also described as the 
exchange of open charm mesons between the charmonium and the incident 
particle. Near the kinematical threshold this exchange was modeled in terms 
of an effective meson Lagrangian . In fig. 1 we show the energy dependence 
of the cross section for J ftp dissociation in this approach. The result for the 
cross section should be considered here as an upper limit since in the model 
the hadronic form-factor was taken to be unity. 

In phenomenological models describing charmonium production in pA 
and light nucleus- nucleus collisions the dissociation cross section of charmo- 
nium on light mesons was assumed to be energy independent 0. Using a 
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Figure 3: As in fig.l. but for the volume of the system: Vhadron is the fraction 
of the volume occupied by the confined-hadronic medium; Vtotai is the total 
volume of the system. 

value of cr^Tv ~ 4.8 mb, which was shown to be consistent with the pA data, 
and the quark structure of hadrons one could fix the o"^(7r,p) ~ 2cr^Ar/3 ~ 3 
mb. 

The differences between these various models for the energy dependence of 
the J/ ip-TT cross section are particularly large near threshold. The theoretical 
uncertainties of the cross section seen in fig. 1 will naturally influence the 
yield of the secondary charmonium production. In the following, we shall 
calculate and compare the yield of the secondary J/tp with all these cross 
sections. 

The ip' absorption cross section has not been evaluated in the quark ex- 
change and effective Lagrangian model. The application of the short distance 
QCD is not possible in this case due to the very low binding energy of the 
ip' ( of the order of 60 MeV) and the correspondingly large size of the ip\ 

For a first estimate of the ip' dissociation on light mesons we assume that 
the absorption cross section is energy independent and attains its geometric 
value of about 10 mb very near the threshold. This assumption was shown 
to provide a quite good description of the ip' suppression observed in S-U 
collisions ||12|| . 
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Figure 4: Time evolution of the number of cc produced in an equilibrium 
QGP for two different values of the initial thermalization time and the initial 
temperature calculated with 12000 pions in the final state. 

2.1 Rate equation for the charmonium production in 
a hadron gas 

In a thermal hadronic medium the rate of charmonium production from 
DD annihilation is determined by the thermally averaged cross section and 
the densities of incoming and outgoing particles. The thermal average of 
the charmonium production cross section < crDfj^^^'^Dj^ > is given by the 
following expression p!3|| : 

where Ki, K2 are modified Bessel functions of the second kind, fn^j^ = 
rriD +mi) and = 771^5 — m^, t = a/s is the center-of-mass energy, /3 the 
inverse temperature, Vat is the relative velocity of incoming mesons and the 
integration limit is taken to be: to = max[{mD + nifj), {m^ + nih)]. 

The thermal cross section for the production and absorption of charmo- 
nium convoluted with the densities of incoming particles describes the rate 
equation for charmonium production in a hadron gas per unit of rapidity by: 
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Figure 5: Time evolution of the number of cc produced in an equilibrium 
QGP for two different values of the charm quark mass. 



where i — {D, D*}, h denotes a tt or p meson and ^ stands for ip' or J/ip. 

In the hadronic medium we include the following secondary charmonium 
production processes: 

D*D + DD* + D*D* ^TT + ip (4) 

D*D* + DD ^ p + ijj (5) 

where ip stands for J/ip or ip'. 

The solution of the rate equation requires additional assumptions on the 
space-time evolution of the hadronic medium and the initial number of D 
and D mesons. 



2.2 Model for expansion dynamics 

We have adopted a model for the expansion dynamics assuming isentropic 
cylindrical expansion of the thermal fireball. To account for transverse ex- 
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Figure 6: Time evolution of the number of cc produced in out of chemi- 
cal equilibrium QGP calculated with the initial conditions from SSPC and 
HIJING. 

pansion we assume a linear increase of the radius of a cylinder with the 
proper time. Thus, at a given proper time r the volume of the system is 
parameterized by: 

V{t) = 7rR^{T)T with R{t) = Ra + 0.15{t - tq) (6) 

where tq is the initial proper time when the system is created as a thermally 
equilibrated quark-gluon plasma and Ra is an initial radius being determined 
by the atomic number of colliding nucleus. For central Pb-Pb collisions 
Ra ~ 6.7 and we fixed tq ~ 0.1 fm. 

The initially produced quark-gluon plasma cools during the expansion 
until it reaches the critical temperature Tc at the time where it starts to 
hadronize. The system stays in the mixed phase until the time r^, where the 
quark gluon plasma is totally converted to a hadron gas. The purely hadronic 
gas can still expand until it reaches the chemical freeze-out temperature Tf 
at the time tj where all inelastic particle scattering ceases and the particle 
production stops. 

The most recent results of Lattice Gauge Theory (LGT) give an upper 



limit for the critical temperature in QCD to be Tc ~ 0.17GeV [0. The 
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Figure 7: Time evolution of Jjil) production from hadron gas due to DD* + 
DD* J /-ifj-K process calculated for Pb-Pb collisions at LHC energy for 
four different parameterizations of the J/ip — it absorption cross sections as 
described in fig. 1. 



analysis of presently available data for different particle multiplicities and 
their ratios measured in heavy ion collisions at SPS energy suggests that the 
chemical freezeout temperature at LHC should be in the range 0.16 < Tf < 



0.18 |T^, [T6|, |T8|, that is very close to Tc. 

To make a quantitative description of the space-time evolution of a ther- 
mal medium one still needs to specify the equation of state. The quark gluon 
plasma is considered as an ideal gas of quarks and gluons whereas the hadron 
gas is described as an ideal gas of hadrons and resonances. We have included 
the contributions of all baryonic and mesonic resonances with a mass of up to 
1.6 GeV into the partition function. To take into account approximately the 
repulsive interactions between hadrons at short distances we apply excluded 
volume corrections. Here we use the thermodynamically consistent model 



proposed in [|T^, |T9| where the thermodynamical observables for extended 
particles are obtained from the formulas for pointlike objects but with the 
shift of the chemical potential. In particular the pressure is described by: 

pextended^j.^ ^ ppointHke ^j. ^ -^^ ^.^j^ /i = - V eigenP'"'"'"^''^ (T , fx) , (7) 



9 



O OQ 




Figure 8: Time evolution of the abundance of if)' from Pb-Pb collisions at 
LHC using constant absorption cross sections for ip' — it and ip' — p oi 10 mb. 

where the particle eigenvolume Veigen = 4|7rr^ simulates repulsive interac- 
tions between hadrons. Following the detailed analysis in [|^, |20[ we have 
assigned the value r = 0.3 fm for all mesons and baryons. 

In fig. 2 we show the time evolution of the temperature for Pb-Pb collisions 
at LHC energy for two different values of T^. The initial temperature of the 
chemically equilibrated plasma was fixed requiring the entropy conservation: 

s,(To)Vo ~ iVj3.6 with = {dNjdy)y=o, (8) 

where A^^ described the final number of pions at midrapidity, Vq is the initial 
volume and Sg(To) the initial entropy density of the thermalized quark-gluon 
plasma. For A^^ = 12000 and ro=0.1 fm we get the initial thermalization 
temperature Tq ~ 1.07GeV for Pb-Pb collisions at LHC. 

The result in fig. 2 shows that for Tc ~ 0.17GeV the QGP starts to 
hadronize after 14 fm and stays in the mixed phase by 20 fm. Taking the 
lowest expected freezeout temperature T/ ~ 0.16 GeV the system exists only 
for a very short time of few fm in a hadronic phase.0 Thus, the contribution 

^we do not count the isentropic expansion phase towards thermal freeze-out, since no 
appreciate particle production takes place there. 
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Figure 9: As in fig. 8 but for different values of the number of DD mesons. 

of the purely hadronic phase to the overall secondary charmonium multiplic- 
ity can be neglected. The secondary charmonium bound states are produced 
almost entirely during the mixed phase. This is in contrast to 0. In fig. 3 we 
show the time evolution of the total volume V as well as the fraction of the 
volume Vh occupied by hadrons during the mixed phase. The total volume 
is calculated following eq.6, whereas Vh can be obtained from the condition 
of entropy conservation in the following form: 

V{t) - V N N 

where and are the entropy density in the quark gluon plasma and the 
hadron gas at the critical temperature. For the equation of state described 
above and T^. = O.lTGeV we have: ~ 11.9 and s^/s^ ~ 0.22. 

The time evolution of the hadronic medium during the mixed phase is 
totally determined by the number of pions in the final state and the entropy 
density of the quark- gluon plasma and the hadron gas at the critical temper- 
ature. Our default parameters for the expansion dynamics are as follows: 

Tc = O.lTGeV , T/ ~ T, , To ~ O.lfm , s" ~ 12/fm^ , sl/sZ ~ 0.22. (10) 
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Figure 10: Total ^' yield produced during the mixed phase as a function of 
pion multiplicity. 

We will study, however, how deviations from the above values influence the 
final multiplicity of the secondary charmonium produced during the evolution 
of hadronic medium. 

2.3 Charm density 

The initial number of open charm mesons at the critical temperature Ndi) 
can be related to the number of primary c, c quarks iVcg produced in A- 
A collisions. Neglecting the possible absorption and production of cc pairs 
during the evolution and hadronization of a quark-gluon plasma we can put 
^Dt) ~ ^cc- Due to charm conservation the number of D's should be equal 
to the number of D mesons. The charm conservation also implies that during 
the mixed phase the number of DD mesons in the hadronic phase N'^j=^is given 
by the fraction of volume occupied by the hadrons that is = NccVh/V 
where V and Vh are described by eq.6 and 9. We further assume that during 
hadronization of the quark gluon plasma the open charm mesons are in a 
local thermal equilibrium with all other hadrons, however, the yield of D 
and D exceeds their chemical equilibrium value. The ratio of the number of 
D to D* mesons at the temperature T is obtained from the relative chemical 
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Figure 11: Total ip' yield produced during the mixed phase as a function of 
the fractions of the critical entropy density of the non-interacting quark-gluon 
plasma. 



equilibrium condition: 



(11) 



The initial number of cc pairs in Pb-Pb collisions at LHC energy was 
obtained by scaling the p-p result with the geometrically allowed number of 
nucleon-nucleon collisions. The cross section for charm production at LHC in 
p-p collisions was calculated by leading order perturbative QCD using code 
PYTHIA with < >= IGeV^ Typical values of the bOcc pairs and 

0.5J/tp were obtained in central Pb-Pb collisions at midrapidity. 

The cc pairs can be also produced and absorbed during the evolution of 
thermally equilibrated quark-gluon plasma. In lowest order in the charm 
quark pairs are produced by gluon and quark pair fusion. Solving a similar 
rate equation as in eq.3 p2| but with the appropriate QCD cross sections we 
calculated in fig. 4 the time evolution of the number of thermal cc pairs. The 
thermal production is seen in fig.4 not to be negligible as compared with the 
preequilibrium production. Dependent on the thermalization time there are 
an additional 5-10 cc pairs produced during the evolution of the thermally and 
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chemically equilibrated plasma. This number, however, depends crucially on 
the value of the charm quark mass rric and the strong coupling constant ag- 
in fig. 5 we compare the thermal cc multiplicity calculated with rric = 1.5 
GeV and rrir = 1.2 GeV with a, = 0.3. 



Dynamical models for parton production and evolution like HIJING 



2^ or SSPC p5[ confirmed, that indeed, after a very short time of the or- 
der of 0.1-0.5fm the partonic medium can reach the thermal equilibrium. 
However, both quarks and gluons my appear far below their chemical satu- 
ration values. Thus, discussing thermal cc production one should take into 
account deviations of gluon and light quark yields from their equilibrium 
values. These deviations can be parameterized by fugacities parameters 
and Xg modifying the Boltzman distribution function of quarks and gluons 

[ii- 

In fig. 6 we calculated the time evolution of cc pairs within the kinetic 



model derived in |23| for a non-equilibrium quark-gluon plasma. The 
following initial conditions for the thermalization time tq, the initial temper- 
ature To and the values for quark Xg and gluon fugacities were fixed for Pb- 



Pb collisions at LHC from SSPC p5[ and HIJING (1) tq = 0.25/m, 

To = 1.02GeV, A° = 0.43 A° = 0.082 (2) nn = 0.5 fm, Tin = 0.82GeV", 
Aj^ = 0.496 A™ = 0.08 H; (3) Tin = 0.5/m, Tin = 0.72GeV, Xf = 0.761 
A*"^ = 0.118 ||2^. In a non-equilibrium plasma, dependent on the initial con- 
ditions, there are 3-5 cc pairs produced during the evolution of plasma as seen 
in fig. 6. One should also note that at the time when temperature reaches 
the critical value of O.lTGeV the quark and gluon fugacities are very close 
to unity which indicates the chemical equilibration of an ideal quark-gluon 
plasma. 

From the results presented in fig's.4-6 one concludes, that dependent on 
the initial conditions as well as on the amount of chemical equilibration of 
the plasma, one expects 5-10 DD pairs at in addition to the 50 DD from 
the hadronization of the initially produced cc pairs. 
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3 Time evolution of the charmonium abun- 
dances 



The number of produced charmonium bound states from DD scattering is 
obtained by multiplying the rate equation eq.3 by the volume of the hadron 
gas eq.9 and then performing the time integration. In fig. 7 we show the time 
evolution of the abundance of J/ip from Pb-Pb collisions at LHC energy 
obtained by solving the rate equation. The calculations were done with the 
four different models for J/ip — ir absorption cross section described in fig.l 
and assuming initially 50 DD and 12000 pions at midrapidity. The results in 
fig. 7 show a very strong sensitivity of J/ip production on the absorption cross 
section. The largest number of the secondary J/ip is obtained with the cross 
section predicted by the quark exchange model. Here the number of J/ip is 
only by 1/2 smaller than the primary value. If, however, the absorption cross 
section is described by short distance QCD the secondary J/ip production is 
lower by more than two orders of magnitude. 

Recent analysis of J/ip photoproduction data confirms the relation be- 
tween the energy dependences of the J/ip cross section and the Feynmann 
x-dependence of the gluon distribution function of the nucleon p, |T0|]. This 
could be used as an indication that the short distance QCD approach is a 
consistent way to calculate the J/ip cross section. Thus, the secondary pro- 
duction of J/ip from the hadronic gas can be entirely neglected in this case. 
The situation can, however, change when discussing the production of ip' 

The time evolution of the abundance of ip' is presented in fig. 8 with the 
same basic parameters as used for J/tp in fig. 7. We show in fig.8 the separate 
contributions of the production processes for the secondary tp' with vr and p 
mesons in the final state. Following the arguments of the previous section the 
absorption cross section for ip' on vr and p mesons were taken to be energy 
independent and equal to its geometric value of 10 mb. The produced ip' 
during the mixed phase from DD annihilation is seen in fig.8 to saturate 
at the large value 0.08, which is almost 1/5 of the primary number of J/ip 
expected for Pb-Pb collisions at LHC energy. 

The yield of the secondary charmonium summarized in fig's. 7-8 depends 
on the parameters used in calculations. In the following we study the modi- 
fication of the results on the secondary ip' by changing the initial number of 
DD mesons as well as the values of the relevant thermal parameters. 



15 



In fig. 9 we discuss the yield of ip' for different multiplicities of DD mesons 
by including the contribution of thermal cc pairs. From the rate equation 
eq.3 it is clear that this dependence is quadratic. The increase of DD by 
10-20% implies the 20 — 40% increase of the ip' yield. 

In the space-time evolution model all parameters were fixed assuming 
= 12000. This number, however, is still not well established and the 
deviations for A^^ in the range 5000 < N^^ < 12000 are not excluded. In fig. 10 
we show the sensitivity of the results for ip' production on A^^r- Keeping the 
same number of DD mesons and decreasing A^^^ the number of the secondary 
i/j' increases as seen in fig. 10. This is mostly because by decreasing A^^r the 
initial density n^j^ of DD increases leading to a larger production of the 
charmonium bound states. 

Modelling the equation of state we have fixed the critical entropy density 
gc-ideai quark gluon plasma by the ideal gas equation of state. From 

the Lattice Gauge Theory we know, however that, due to interactions, the 
entropy density could be smaller by a significant factor. Deviation of the 
momentum distribution of quarks and gluons from the chemical saturation is 
also reducing the critical entropy density of quark-gluon plasma. To estabhsh 
the influence of a decreasing critical entropy density we calculated the total 
number of ip' as a function of the ratio s^/s^"*'^'^"' in fig.ll. Decreasing the 
entropy density of a quark-gluon plasma by a factor of two reduces the yield 
of J/ip by 70%. This is mostly because the time the system spends in the 
mixed phase is shorter. However, even with this reduction of the entropy 
density the secondary production of ip' is not negligible. The influence of the 
result on the critical entropy density in a hadron gas is less important. 
Increasing s^, by a factor of two reduces the total number of ip' by only 10%. 

The secondary ip' production in a mixed phase is also sensitive to the 
parameterization of the time evolution of the system size in transverse di- 
rection. Increasing the transverse expansion parameter from 0.15 to 0.45 in 
eq.6 decreases the total number of ip' produced during the mixed phase by 
40 %. Taking a quadratic dependence of -R(r) on proper time as proposed 
in [0] decreases the yield of ip' from 0.08 to 0.036, which is still comparable 
with the number of primary ip'. 
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4 Conclusions 



We have considered the possibihty of the secondary charmonium production 
in ultrarelativistic heavy ion colhsions at LHC energy. Admitting thermal- 
ization of a partonic medium crated in a colhsion and the subsequent first 
order phase transition to a hadronic matter we have shown that the sec- 
ondary charmonium production appears almost entirely during the mixed 
phase. The yield of secondarily produced ip mesons is very sensitive to the 
hadronic absorption cross section. Within the context of the short distance 
QCD approach this leads to negligible values for J/ip regeneration. The 
production, however, can be large and may even exceed the initial yield from 
primary hard scattering. Thus it is conceivable that at LHC energy the ■0' 
charmonium state can be seen in the final state whereas J/ production can 
be entirely suppressed. The appearance of the -0' in the final state could be 
thus considered as an indication for the charmonium production from the 
secondary hadronic rescattering. 
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